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weighted  forcing  function  for  i^  mode 

radius  of  circle 

time 

sampling  time 

frequency  response  function  of  pitch  with  respect 
to  wave 

horizontal  fluid  velocity 
vertical  fluid  velocity 
ship  forward  speed 
shear  force 

vertical  wave  orbital  velocity 
relative  vertical  velocity 

horizontal  axis  in  direction  of  forward  motion 
of  ship  (along  length  of  ship) 

mode  shape  of  i*"*1  mode 

horizontal  axis  directed  to  starboard 

heave  motion,  positive  upwards 

vertical  elastic  deflection 

relative  vertical  immersion  change 

elastic  deformation  angle 

surface  wave  elevation,  positive  upwards 

surface  wave  elevation  encountered  while  moving  forward 

pitch  angle,  positive  bow  up 

wavelength 
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sum  of  sectional  ship  mass  and  added  mass 
weighted  total  mass  for  ifc^  mode 
velocity  potential 
velocity  potential  of  surface  waves 
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INTRODUCTION 


When  considering  the  structural  loads  and  responses  acting 
on  a  ship  in  a  seaway,  the  different  effects  contributing  to  the 
midship  vertical  bending  moment  (which  is  the  major  hull  girder 
structural  load)  must  be  recognized  and  treated  separately  in  any 
analysis  and/or  design  study.  In  the  case  of  the  vertical  bend¬ 
ing  moment  on  a  ship  in  waves,  the  total  load  is  made  up  of  two 
distinct  contributions,  viz.  the  slowly- varying  bending  moment 
directly  induced  by  the  waves,  and  also  a  higher  frequency  response 
that  exhibits  the  vibratory  characteristics  associated  with  the 
structural  modes  of  the  ship.  The  slowly-varying  bending  moment 
has  a  frequency  domain  bandwidth  that  lies  within  that  associated 
with  the  waves  (i.e.  as  encountered  by  the  ship  in  its  forward 
motion  at  speed) ,  while  the  vibratory  response  is  of  a  much  higher 
frequency  that  is  most  often  associated  with  the  first  structural 
bending  mode  (2-noded  vibration)  or  perhaps  some  of  the  higher 
modes  (see  [1]). 

Methods  of  computing  the  wave-induced  bending  moments,  both 
vertical  and  lateral,  as  well  as  the  wave-induced  torsional  moment, 
have  been  developed  within  the  course  of  a  project  sponsored  by 
the  Ship  Structure  Committee  (Project  SR-174,  "Ship  Computer 
Response"),  and  the  description  of  the  analytical  methods,  cal¬ 
culation  results,  comparison  with  model  experiments,  and  the  developed 
digital  computer  programs  are  provided  in  [2]  and  [3]  .  These  responses 
have  been  found  to  be  linear  in  regard  to  their  variation  with  wave 
amplitude;  they  are  continuously  varying  just  as  the  waves  vary  con¬ 
tinuously;  their  statistical  properties  are  determined  via  spectral 
analysis  techniques,  with  these  properties  being  statistically 
stationary  in  the  same  manner  as  the  waves;  and  their  evaluation 
requires  determining  the  rigid  body  responses  of  the  ship  hull  to¬ 
gether  with  the  contribution  of  these  rigid  body  motions  to  the  local 
forces  that  are  distributed  along  the  ship  hull  (outer  envelope 
within  the  water  as  well  as  the  internal  ship  mass  distribution). 

The  developments  in  [2]  and  [3]  provide  a  useful  tool  for  valid 
estimation  of  bending  moment  responses  directly  induced  by  the  waves, 
with  the  output  in  terms  of  statistical  averages  appropriate  to  the 
particular  wave  spectrum  in  which  the  ship  under  consideration  is 
assumed  to  operate. 

The  vibratory  structural  responses,  as  reflected  in  the 
vertical  bending  moment,  are  caused  by  entirely  different  mechanisms, 
and  the  properties  of  such  responses  also  differ  from  those  of  the 
slowly  varying  wave-induced  bending  moments.  These  responses  are 
often  associated  with  the  occurrence  of  large  ship  motions  where 
emersion  of  the  bow  region  can  occur,  leading  to  impact  forces 
associated  with  entering  the  water  (i.e.  the  ordinary  case  of  ship 
slamming  phenomena)  or  in  the  case  of  other  ships  with  large  bow 
flare  present  that  leads  to  "whipping"  resulting  from  the  forces 
developed  due  to  the  bow  flare  shape  variations  (e.g.  see  [1]).  The 
force  mechanisms  associated  with  the  ship  motions  that  produce  the 
resulting  input  excitation  are  thus  dependent  on  nonlinear  effects. 


although  the  rigid  body  motions  of  the  ship  per  se  are  suf¬ 
ficiently  well  represented  as  linear  responses,  and  these  forces 
are  known  to  be  impulsive  in  nature.  The  resulting  structural 
response,  due  to  the  excitation  of  the  basic  structural  modes 
of  the  ship,  is  then  manifested  as  a  series  of  non-continuous  high 
frequency  oscillations  (i.e.  in  the  vertical  bending  moment)  that 
only  occur  following  the  development  of  the  impulsive  local  forces 
at  the  bow,  i.e.  a  nonstationary  record.  The  frequency  of  these 
vibratory  responses  is  usually  that  of  the  first  structural  mode, 
i.e.  the  two-noded  vibratory  mode  {see  [1]),  and  the  oscillations 
decay  as  a  result  of  the  combined  influence  of  both  structural 
and  hydrodynamic  deunping.  Thus  it  is  easily  seen  how  these  struc  - 
tural  responses  due  to  sleunming  phenomena  differ  substantially  from 
those  bending  moments  that  are  the  direct  wave-induced  structural 
reactions . 

In  addition  to  the  vibratory  bending  responses  that  arise 
from  sleunming  effects  (both  bow-flare  as  well  as  bottom  impact 
sleunming) ,  another  source  of  such  responses  occurs  when  the  ship 
has  small  (or  insignificant)  motions  in  relatively  short  waves  such 
that  the  frequency  of  encounter  with  the  waves  is  close  to  that  of 
the  lower  structural  modes  of  vibration  of  the  ship.  This  particu¬ 
lar  phenomenon  has  been  denoted  as  "springing",  with  recent  analyti¬ 
cal  efforts  aimed  at  studying  this  effect  given  in  [4]  and  [5]  .  The 
particular  ships  for  which  springing  has  an  important  influence  are 
large  tankers  and  bulk  carriers  (such  as  illustrated  in  [6]),  as 
well  as  for  fast  ships  such  as  destroyers  and  container  ships. 
According  to  the  results  in  [4]  and  (5)  ,  as  well  as  the  physical 
interpretation  of  the  influence  of  short  waves,  i.e.  it  is  a  direct 
wave-excited  vibration  that  is  dominated  by  the  two-node  vibratory 
bending  resnonse  of  the  ship  structure.  The  application  of  spectral 
analysis  techniques  to  this  linear  response,  as  shown  in  [4]  and  [5] 
provides  statistical  measures  of  the  bending  moment  due  to  springing 
These  values  can  then  be  compared  to  the  values  due  to  the  direct 
wave-induced  bending  moments  experienced  under  the  same  and/or 
other  particular  wave  spectral  conditions  in  order  to  assess  magni¬ 
tudes  that  would  be  useful  for  design  purposes. 

It  must  be  recognized  however,  that  the  springing  phenomenon 
involves  consideration  of  short  waves  relative  to  the  ship  length, 
and  that  the  theoretical  bases  for  evaluating  the  wave  excitation 
forces  acting  on  a  ship  were  originally  derived  for  conditions  where 
the  wavelengths  were  long  compared  to  the  'imensions  of  the  ship 
cross-section  (i.e.  beam  and/or  draft).  As  a  result  there  remains 
some  question  as  to  the  full  validity  of  any  results  obtained  from 
such  analyses,  at  least  until  some  further  consideration  has  been 
given  to  a  more  precise  analysis  of  the  forces  acting  on  a  ship  hull 
due  to  short  waves. 

In  order  to  obtain  information  on  the  vibratory  structural 
loadings  associated  with  the  various  phenomena  described  above, 
the  techniques  of  computer  simulation  can  be  extended  to  these 
cases,  based  upon  some  of  the  developments  described  previously  in 
[1J.  In  view  of  the  fact  that  the  slamming  and  bow  flare  effects 
require  nonlinear  *orce  determination,  as  well  as  the  fact  that 
the  resu1  .ing  bending  moment  variations  are  also  nonstationary,  a 


time  domain  treatment  is  suggested  (see  [1]  and  [7]).  Similar 
treatment  in  regard  to  a  time  domain  output  representation  can 
also  be  applied  to  the  case  of  springing,  although  it  is  possible 
to  obtain  data  from  the  use  of  frequency  response  and  spectral 
analysis  methods  for  that  particular  effect.  However  the  form  of 
time  domain  outputs,  comparing  the  slowly-varying  wave-induced 
vertical  bending  moment  with  the  vibratory  bending  moment,  are  al¬ 
ways  useful  illustrations  that  can  provide  certain  additional  in¬ 
sight.  Another  possible  use  of  time  domain  outputs  would  be  a 
more  direct  method  of  comparison  with  experimental  results,  since 
that  would  be  a  very  definite  method  of  validation  of  any  analyti¬ 
cal  procedure,  as  long  as  the  complete  input  information  required 
for  such  a  comparison  is  available. 

Since  the  presence  of  time  histories  for  slamming-related 
phenomena  is  the  only  possible  way  of  providing  a  proper  analysis 
of  these  effects,  the  ability  to  generate  time  histories  with  a 
computer  at  a  relatively  fast  rate  as  compared  to  real  time  (i.e. 
the  actual  time  for  such  occurrences  as  recorded  on  full  scale 
ships)  is  another  useful  characteristic  of  computer  simulation  as 
compared  to  obtaining  full  scale  data  and/or  carrying  out  model 
tests  with  proper  structural  simulation  in  the  model.  Thus  the 
proper  development  of  computer  simulation  for  these  vibratory 
structural  loadings,  which  are  assuming  more  importance  as  ships 
become  larger  and  faster,  will  provide  a  useful  tool  in  design, 
analysis,  and  evaluation  studies  associated  with  modern  ship  develop¬ 
ment. 

The  present  report  is  intended  to  describe  the  results  of 
computer  simulation  techniques  for  determining  the  various  types 
of  vibratory  structural  loads  experienced  by  a  ship  at  sea,  as 
represented  by  the  vertical  bending  moment.  A  number  of  particular 
computations  will  be  illustrated  and  the  nature  of  the  results  com¬ 
pared  with  that  experienced  in  model  tests  and/or  full  scale  tests, 
in  accordance  with  the  extent  of  available  data  for  such  purposes. 

An  evaluation  of  the  computational  procedures,  with  regard  to  time 
and  cost,  will  also  be  provided  in  order  to  illustrate  the  possible 
benefits  that  can  be  obtained  via  computer  simulation.  These  re¬ 
sults  do  nbfc  represent  the  final  capability  and  results  of  computer 
simulation,  but  will  only  serve  to  illustrate  the  prospect  of  apply¬ 
ing  such  techniques  in  a  more  comprehensive  manner  that  includes  a 
more  complete  mathematical  representation  and  that  would  provide 
a  final  computer  program  output  that  could  be  applied  toward  the 
routine  evaluation  of  these  particular  structural  loadings,  in  the 
same  manner  as  was  provided  in  [2]  and  [3], 

MODAL  METHOD  OF  ANALYSIS 

The  vibratory  structural  loads  are  determined  by  considering 
the  ship  structure  to  be  an  elastic  beam  with  nonuniform  mass  and 
elastic  beam  with  its  length.  The  equations  representing  this 
type  of  model  have  been  presented  previously  in  [1] ,  and  they  are 
given  below  as: 

a7z„  3  z  3  V 

u  - -  ♦  c  — £  ♦  — 2-  «  P(x,t)  ,  (1) 

at2  at  3x 
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where  m  *=  n  (x)  is  the  sum  of  the  ship  mass  and  the  added  mass  at 
a  section;  z  represents  the  vertical  elastic  deflection;  c  is 
the  damping  coefficient;  V  is  the  shear  force;  and  P(x,t)  is  the 
local  input  force  due  to  snip-wave  interaction. 


—  «  V_  +  I 


where  M  is  the  bending  moment,  I  is  the  mass  moment  of  inertia  of 
a  section,  and  y  is  a  deformation  angle,  with  the  last  term  on  the 
right  in  Equation  (2)  representing  the  rotary  inertia. 


M  =  El 


is  the  fundamental  elastic  equation,  with  El  the  bending  flexural 
rigidity. 


relates  the  bending  and  shear  effects,  where  KAG  is  the  vertical 
shear  rigidity. 

These  partial  differential  equations  were  considered  in  [1] , 
where  different  possible  procedures  for  solution  were  examined  for 
the  case  of  an  impulsive  force  localized  in  Uie  bow  region  (the 
problem  of  bow-flare  slamming) .  Since  the  output  in  time  history 
form  was  desired,  for  a  relatively  higher  frequency  phenomenon 
relative  to  the  direct  wave  system  and  its  slowly-varying  ship 
responses,  the  various  methods  of  direct  solution  of  the  partial 
differential  equations  by  conversion  into  ordinary  differential- 
difference  equations,  breaking  up  the  beam  into  a  large  number  of 
nodal  segments,  etc.  were  found  to  be  inapplicable  to  the  simulation 
requirements.  This  was  based  upon  limitation?  inherent  in  the 
computer  components,  problems  of  computational  "stability"  of  the 
solutions,  as  well  as  the  basic  desire  to  obtain  solutions  at  rates 
significantly  faster  than  real-time  (see  [1]  for  a  detailed  dis¬ 
cussion  of  these  different  aspects  of  computer  simulation) . 

The  method  that  was  applied  in  [1]  for  treating  the  general 
problem  of  vibratory  structural  responses  was  to  use  a  modal  model 
to  represent  the  basic  ship  beam  structure,  with  neglect  o7“rotary 
inertia  (which  should  have  a  negligible  effect  for  the  present  class 
of  applications) .  The  different  variables  in  the  equations  are 
represented  in  product  form  as 


*e(x,t) 
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M  (x,  t) 


00 

-I* 


(t)Mi (x) 


5 


(5) 


vs(x't)  ■  L  vt)v.  <x> 

i-1  1 

where  X.  (x)  is  the  normal  mode  shape  of  the  mode.  Using  the 
resultsxof  a  separation  of  variables  method  of  solution  for  the 
unforced  beam  motion  solution  as  a  basis,  the  forced  motion  responses 
were  found  to  be  represented  by 

Wi^  +  +  Kiqi  «  Qi(t)  (6) 


where  L/2 

vi  =*  wXi2(x)  dx 

•  -L/2 


(7) 


C.  -  v. 


wi  \V\  Mi 

since  c/y  is  assumed  to  be  constant  along  the  ship  (at  least  for 
the  structural  damping  contribution)  , 


Ki  *  wi2wi  ' 

±.y. 

with  u.  the  natural  frequency  of  the  i  mode  (in  rad./sec.),  and 
1  (L/2 


(t) 


P(x,t)Xi(x)  dx 


'-L/2 


The  representation  of  the  bendi.ig  moment  spacial  weight¬ 
ing  is  given  by 


Mi  "  wi2 


(x-s)y (s)X^ (s)ds 


-L/2 


where  the  position  x  is  the  location  at  which  the  bending  moment 
is  desired,  with  x»0  the  midship  location  taken  as  the  coordinate 
origin  and  all  integrations  over  the  ship  length  extending  from 
x«  -L/2  (stern)  tc  x  -  L/2  (bow).  The  expression  in  Equation  (11) 
is  based  upon  considering  the  main  contribution  to  the  bending 
moment  to  arise  from  the  resulting  inertial  loads  along  the  hull 
due  to  the  vibratory  deflections,  including  the  fluid  inertial  force 
associated  with  the  added  mass. 


The  solution  of  Equation  (6)  is  to  be  obtained  for  each 
mode,  and  then  weighted  in  accordance  with  the  particular  mode 
shape  or  related  spacial  function,  as  indicated  in  Equation  (5) 
where  all  responses  are  represented  as  the  sum  of  the  individual 
responses  excited  in  the  different  modes  of  vibration.  However 
in  actual  practice,  for  the  ship  responses  considered  herein,  the 
first  mode  of  vibration  is  most  predominant  with  only  a  negligible 
contribution  from  the  higher  modes,  thereby  simplifying  the 
computational  requirements .  In  order  co  carry  out  the  computer 
simulation  by  this  method  it  is  necessary  to  establish  procedures 
for  determining  the  natural  frequencies  and  mode  shapes  (eigenvalues 
and  eigenfunctions)  associated  with  a  ship  structure,  and  also  the 
technique  for  representing  and/or  determining  the  local  input  forces 
due  to  ship-wave  interaction. 

The  determination  of  the  natural  frequencies  and  mode  shapes 
of  a  ship  structure,  represented  as  a  free-free  beam,  is  carried 
out  by  adapting  the  available  results  of  the  U.  S.  Navy  Generalized 
Bending  Response  Code  (GBRC1) ,  as  described  in  [8].  When  input 
data  in  the  form  of  the  distribution  of  inertial  and  structural 
properties  (such  as  bending  stiffness,  shear  stiffness,  etc.)  is 
given,  a  digital  computer  program  based  on  Equations  ( 1 )  —  ( 4 )  will 
provide  the  desired  output.  The  digital  computer  program  described 
in  [8]  was  obtained  and  converted  for  use  on  the  CDC  6600  digital 
computer,  and  a  check  on  the  results  was  obtained  by  comparing 
the  output  for  the  first  mode  frequency  and  mode  shape  of  the  USS 
ESSEX  with  that  given  in  [9] ,  which  was  used  in  the  work  of  [1] . 

Very  good  agreement  was  obtained  in  that  case,  and  hence  a  useful 
tool  for  determining  ship  natural  frequencies  and  mode  shapes  is 
thus  available  as  one  element  for  use  in  treating  problems  of  ship 
structural  response  in  waves. 

The  evaluation  of  the  external  forcing  function  is  another 
problem,  and  that  depends  upon  the  basic  cause  of  the  particular 
type  of  vibratory  loading  under  consideration.  In  the  case  of  bow 
flare  slamming  treated  in  11],  the  force  considered  as  the  input 
to  excite  vibratory  response  was  due  to  nonlinear  variations  in 
buoyancy  and  inertial  forces,  over  and  above  those  used  in  the 
linear  ship  motion  analysis.  The  instantaneous  immersion  and 
relative  velocity  of  the  various  ship  sections  determines  the 
buoyancy  and  fluid  momentum  values,  using  tabulated  values  of 
sectional  area  and  added  mass  obtained  from  a  previous  off-line 
computation  (see  [1]).  In  the  case  of  slamming  that  involves 
bottom  impact  when  bow  sections  re-enter  the  water  after  emer¬ 
sion,  the  force  computation  also  involves  determination  of  non¬ 
linear  variations  of  added  buoyancy  and  inertial  contributions 
that  are  related  to  instantaneous  values  of  section  immersion 
and  relative  velocity,  acceleration,  etc.  As  a  result,  in  both 
of  these  cases,  it  is  necessary  to  obtain  time  history  representa¬ 
tions  of  the  different  linear  ship  motions  necessary  for  evalua¬ 
ting  the  different  force  terms  described  above. 

This  is  accomplished  by  means  of  a  convolution  integral 
operation  in  the  time  domain,  where  the  encountered  wave  motion 
time  history  at  a  reference  point  relative  to  the  ship  is  the 
input  data  that  is  operated  upon  with  a  weighting  function  kernel. 
The  kernel  function  is  obtained  as  the  Fourier  transform  of  the 
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frequency  response  function  for  the  variable  of  interest,  e.g. 
in  the  case  of  pitch  motion  the  time  history  is  given  by 


e(t) 


K  (t-x) «  (t)<3t 
9  m 


(12) 


—  ao 


where  n  (t)  is  the  encountered  wave  motion  time  history  (present 
and  pas?  record) , 


Ke(t) 


00 


Ten(“e;xl) 


(13) 


is  the  pitch  kernel  function,  and  the  pitch  frequency  response  to 
a  unit  sinusoidal  wave  (as  measured  at  a  point  x.  ahead  of  the 
origin  of  coordinates  on  the  ship)  is  given  by 


Te„<Vxi>' 


M 
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sl[vF(ue)xl 


(14) 


in  terms  of  amplitude  (relative  to  wave  amplitude  a)  and  phase 
Ue,  relative  to  the_coordinate  origin,  with  the  correction  to  x. 
in  terms  of  F(u  )  *  — j  represented  as  a  function  of  a>  for  a  fixed 
forward  speed) .  These  procedures  are  described  and  illustrated  in 
[1]  for  a  particular  ship  case,  and  other  illustrations  of  this 
method  as  well  as  extensions  to  determine  further  aspects  of  time 
domain  computations  are  given  in  [10] . 

All  of  these  operations  are  based  on  application  of  digital 
computers  to  carry  out  the  evaluation  of  the  frequency  responses , 
kernel  functions,  etc.,  as  well  as  the  evaluation  of  the  convolution 
integral  operation  on  a  wave  motion  time  history  input  to  produce 
outpr„  time  histories  of  the  desired  motion,  etc.  The  digital 
computer  would  then  provide  the  necessary  hydrodynamic  force  time 
history,  serving  as  a  function  generator  essentially,  in  a  hybrid 
computer  simulation  procedure  where  the  analog  computer  portion 
would  provide  the  solution  of  the  elastic  response  equation  in 
modal  form  given  by  Equation  (6).  The  required  hybrid  linkage 
elements  such  as  A-D  converters,  D-A  converters,  etc.,  together 
with  a  noise  generator  and  wave  spectrum  shaping  filter,  when  added 
to  the  digital  and  analog  computers  described  above,  form  the 
complete  hybrid  computer  simulation  system  that  was  proposed  and 
applied  in  [1]  with  good  success. 

Considering  the  greater  generality  and  availability  of  large 
digital  computers,  as  well  as  the  fact  that  the  slowly-varying  wave- 
induced  bending  moments  as  well  as  frequency  response  functions  of 
ship  motions  required  for  vibratory  responses  are  found  by  use  of 
large  digital  computers,  the  possibility  of  solving  for  the  vibratory 


structural  responses  in  time  history  form  on  such  computers  would  be 
a  useful  simplification  of  computational  requirements.  One  possible 
approach  would  be  to  evaluate  the  closed  form  solution  of  Equation 
(6) ,  which  would  be  expressed  as 


q^t) 


where 


Qi(x)  -fe(t-T) 
-  e  si 


o  Vi 


sin  >^(t-i)dT 
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which  is  the  solution  for  the  initial  conditions  q . (0)=q. (0)=0. 
However  a  much  simpler  method  of  solution  is  available  wnen  re¬ 
cognizing  that  the  differential  equation  in  Equation  (6)  can  be 
represented  as  a  recursive  digital  filter,  as  illustrated  in  [11], 
where  the  term  "recursive"  implies  that  the  computation  of  the 
output  is  an  explicit  function  of  previous  outputs  and  inputs.  The 
general  second  order  difference  equation  for  this  system  is 


q (nT) =A^q (nT -T)+A2q (nT-2T)+G (nT) 


where  T  is  the  sampling  period  and  the  coefficients  A^  and  A^  are 
related  to  the  coefficients  in  Equation  (6) . 

In  order  to  check  the  capability  of  this  digital  model  to 
represent  the  required  solutions,  the  case  with  coefficients 
corresponding  to  the  first  mode  of  the  USS  ESSEX  treated  in  [1] 
was  established  with  an  input  function  represented  by  the  output 
of  bandpass  filtered  white  noise  signal.  A  pulse  was  applied  to 
this  noise  signal  (which  contained  frequencies  substantially 
lower  them  the  representative  second  order  dynamic  system) ,  and 
the  response  was  obtained  simultaneously  from  an  analog  computer 
and  from  a  digital  computer  operating  with  the  required  A-D  and  D-A 
converters  necessary  for  computing  and  display  purposes.  Repre¬ 
sentative  output  signals  from  both  computers,  as  well  as  the  input 
signal  used,  are  shown  in  Figure  1  where  it  can  be  seen  that  the 
digital  simulation  is  essentially  the  same  as  the  results  of  analog 
simulation.  This  simple  experiment  was  performed  on  a  small 
digital  computer  (PDP-8)  using  a  sampling  rate  of  12/sec.,  thereby 
indicating  an  expectation  of  excellent  digital  simulation  by  this 
technique  when  using  a  larger  (and  faster,  more  accurate,  etc.' 
digital  computer. 

The  application  of  the  methods  described  in  this  section  to 
the  different  types  of  vibratory  structural  responses  of  a  ship 
in  waves  is  presented  in  the  following  sections  of  this  report. 

The  particular  force  mechanism  proposed,  the  method  of  time-domain 
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representation,  and  the  results  obtained  are  described  separately 
for  each  type  of  response  phenomenon.  All  necessary  simulation 
techniques  used  are  described  together  with  information  on  the 
time  requirements  as  compared  to  the  extent  of  real  time  simulated. 


•log  computer  output 


digital  ooaputar  output 


\  noiao  ♦  pulaa  signal  input 

Fig.  1.  Comparison  of  Output  Sig¬ 
nals  from  Analog  and  Dig¬ 
ital  Computer  Simulations 


APPLICATION  TO  BOW  FLARE  SLAMMING 


The  case  of  bow  flare  slamming  has  previously  been  treated 
in  [1],  and  the  present  set  of  computer  experiments  is  being  carried 
out  in  order  to  judge  the  capabilities  of  a  pure  digital  simulation 
of  the  bending  moment  due  to  the  effects  of  large  bow  flare,  where 
the  treatment  here  is  restricted  to  the  case  of  head  seas.  The 
input  force  is  made  up  of  two  terms,  an  inertial  term  represented  by 


Px(x,t)  * 


Dt  (ranlwr) 


where  the  operator 
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with  V  the  ship  forward  speed,  m  is  the  additional  added  mass 
at  a  section  that  is  determined  fPfim  the  instantaneous  immersion 


•  (7 


geometry  of  the  ship  section  after  subtracting  out  the  added  mass 
determined  from  the  still  water  (linear  theory)  reference  geometry, 
and  wr  is  the  relative  velocity  at  the  section,  given  by 

wr  =  z  +  xe  -  ve  -  WQ(x,t)  (19) 


The  rigid  body  motions  z  and  e  (and  their  derivatives)  are 
determined  from  linear  theory  solutions  (see  [2]),  and  w  (x,t) 
is  the  wave  orbital  velocity  given  by  ° 

wQ(x,t)  =  ^  cos  'p-  +  (V+c)t]  (20) 


for  'he  present  head  sea  case  (illustrated  here  for  sinusoidal 
waves),  where  c  is  the  wave  propagation  speed.  The  force  due 
to  buoyancy,  denoted  as  P2(x,t),  is  represented  by 

P2(x,t)  =  pgAn£(zr;x)  (21) 


where  K  is  the  additional  cross-sectional  area  at  a  section  due 
to  the  Sifference  between  the  area  corresponding  to  the  instan¬ 
taneous  submerged  portion  of  the  ship  section  and  that  corresponding 
to  the  still  waterline,  after  eliminating  the  linear  buoyancy  force 
terms.  The  quantity  A  is  determined,  for  a  particular  ship 
section,  as  a  functionnof  the  relative  immersion  change 


z  =  z  +  x0  -  n (x,t)  (22) 

• 

and  it  is  expressed  as 


V  -  A  -  Ao  +  B*zr  (23> 

where  A  is  the  instantaneous  submerged  area  of  a  section,  A  is 
the  area  up  to  the  still  waterline,  and  B*z  corresponds  to°the 
linear  spring  rate  that  is  included  in  the  determination  of  the 
direct  wave-induced  rigid  body  motions  and  the  wave-induced  verti¬ 
cal  bending  moment  (B*  is  local  beam) . 

The  nonlinear  buoyancy  force  defined  in  Equations  (21) -(23) 
is  determined  in  tabular  form  at  various  stations  from  the  ship 
lines  drawing,  and  the  values  of  the  added  mass  for  different  ship 
sections  are  calculated  for  the  different  levels  of  immersion.  The 
added  mass  used  in  this  investigation  is  the  high  frequency  limit 
appropriate  to  vibratory  response,  and  it  is  independent  of  gravity 
wave  effects  and  is  hence  frequency-independent.  A  generalized 
program  for  computing  the  two-dimensional  added  mass  of  arbitrary 
ship  sections  has  been  developed  and  described  in  [1). 


The  computations  illustrating  bow  flare  slamming  are  made 
for  the  USS  ESSEX,  using  the  previously  determined  frequency  re¬ 
sponse  characteristics  and  resulting  kernel  functions  for  relative 
immersion  and  relative  immersion  velocity  obtained  in  the  work  of 
[1] .  These  quantities  were  determined  for  a  13.8  knot  forward 
speed  and  the  reference  position  at  which  the  waves  were  measured 
(the  location  of  x.)  was  taken  at  30  ft.  ahead  of  the  ship  FP.  This 
location  was  sufficiently  far  forward  so  that  the  kernel  functions 
had  no  significant  magnitude  for  negative  values  of  their  argument, 
which  is  necessary  for  producing  a  kernel  function  that  would  allow 
evaluation  of  instantaneous  conditions  without  any  time  lags,  i.e. 
the  operation  takes  place  only  on  present  and  past  values  of  the 
wave  motion  time  history.  Illustrations  of  the  frequency  responses 
and  time  domain  kernel  functions  for  this  case  are  given  in  [1], 
where  a  time  domain  representation  is  also  available  in  the  form 
of  a  convolution  integral  for  the  slowly- varying  bending  moment 
directly  induced  by  the  waves.  This  quantity  i*  represented  by 


M^t)  =  KM(t)nm(t-T)dT  (24) 


where  K„(t)  is  found  as  the  Fourier  transform  of  the  bending  moment 
frequency  response  (i.e.  amplitude  and  phase). 

The  wave  record  that  is  used  in  the  digital  simulation  is 
obtained  by  initially  constructing  a  digital  filter  that  would 
produce  a  fit  to  the  form  of  power  spectral  density  functions  that 
represent  the  wave  spectrum  as  encountered  by  the  ship  when  moving 
forward  at  speed,  i.e.  a  spectrum  in  terms  of  the  frequency  of 
encounter  w  .  A  sequence  of  random  numbers  that  represent  the 
equivalent  of  white  noise  is  passed  through  this  filter  to  produce 
the  desired  representation  of  the  wave  spe strum,  in  accordance  with 
the  procedures  described  in  [12]  to  produce  pseudo-random  sequences 
with  limited  bandwidth  that  represents  a  simulation  of  the  wave 
random  process.  This  method  of  fit  to  the  wave  spectrum  matches  the 
general  form,  the  location  of  the  frequercy  of  maximum  spectral 
power,  and  the  rms  value  (or  significant  height)  of  the  resulting 
time  history. 

The  nonlinear  hydrodynamic  forces  were  determined  and  applied 
at  four  points  of  the  twenty  beam  segments  making  up  the  USS  ESSEX. 
They  were  applied  at  the  midpoints  of  each  of  the  first  four  seg¬ 
ments  (stations  19 18 J,  17 j  and  16 J)  with  the  appropriate  weight¬ 
ing  in  terms  of  the  mode  shape  according  to  Equation  (10).  Computer 
runs  were  made  to  simulate  the  bending  moment  responses  using  a 
CDC  6600  digital  computer,  where  the  runs  were  made  for  a  very  severe 
sea  3tate,  viz.  Sea  State  9  which  corresponded  to  a  wind  speed  of 
50  knots.  Separate  outputs  were  obtained  on  the  computer  for  the 
midship  bending  moments  due  to  waves  and  due  to  bow  flare  slamming 
(for  the  bending  moment  at  midship) .  The  output  for  the  bow  flare 
slam-induced  bending  moment  was  taken  to  be  that  due  to  the  first 
mode  response  only,  and  all  other  higher  modes  were  neglected  in 
accordance  with  previous  results  indicating  validity  of  this  proce¬ 
dure  in  [1] . 
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The  output  in  the  form  of  time  histories  of  the  wave-induced 
midship  bending  moment  and  the  total  bending  moment  (sum  of  wave- 
induced  and  bow  flare  slam-induced  bending  moments)  is  shown  in 
Figure  2  together  with  the  wave  record  corresponding  to  this  Sea 
State  9.  A  relatively  short  time  interval  of  only  17  seconds  is 
illustrated  in  Figure  2  since  the  computations  were  carried  out  in 
order  to  correspond  to  every  0.5  seconds  of  real  time,  and  extensive 
plotting  is  required  for  a  relatively  short  time  extent.  However, 
the  figure  amply  illustrates  the  presence  of  the  large  slamming 
effect  evidenced  by  the  bending  moment  response  due  to  the  "whipping" 
associated  with  the  large  bow  flare  slamming.  While  this  is  a  very 
severe  case,  the  capability  of  representing  this  response  with  a 
digital  computer  is  amply  illustrated  by  these  results.  The  magnitude 
of  the  maximum  total  bending  moment  double  amplitude  at  midship, 
corresponding  to  the  78  ft.  significant  height  of  the  waves,  was  found 
to  be  5  x  10**  ft. -tons,  which  is  larger  in  this  case  than  the  values 
obtained  in  the  earlier  work  in  [1]  .  This  was  due  to  a  more  precise 
inclusion  of  the  effects  of  the  spacial  derivative  term  in  the  inertial 
force  defined  in  Equation  (18),  i.e.  the  term  corresponding  to  . 
which  produces  results  that  are  closer  to  available  model  test  3x 
data  for  this  condition.  The  good  agreement  in  overall  character¬ 
istics  of  response,  as  well  as  the  magnitude  obtained  in  [1]  and 
the  present  simulation  results  as  compared  to  model  test  data,  shows 
that  the  simulation  provides  generally  good  prediction  of  expected 
structural  responses  of  ships  in  waves  when  experiencing  bow  flare 
slamming. 
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Fig.  2.  Midship  Bending  Moment  Time  Histories,  USS  ESSEX, 

V  «  13.8. kts..  Sea  State  9 

In  carrying  out  the  various  computational  steps  required  to 
evaluate  the  bending  moment  time  histories  for  particular  sample 
sea  states,  some  of  the  associated  computations  can  be  considered 
separately  from  the  actual  evaluation  of  the  total  bending  moment 
(sum  of  wave-induced  and  bow  flare  slam-induced  effects)  per  se.  Thus 
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'■-he  computation  of  the  structural  mode  shape  and  frequency,  the 
kernel  functions  for  relative  immersion  and  velocity,  the  deter¬ 
mination  of  the  tabulated  values  ol  buoyancy  and  added  mass,  the  wave 
spectrum  filter,  etc.  can  be  carried  out  just  once  for  any  particular 
ship  at  a  given  forward  speed  since  these  quantities  are  then  uni¬ 
versal  for  those  particular  conditions.  Once  they  are  established, 
the  actual  time  history  of  bending  moments  is  then  determined  by 
the  procedures  discussed  previously  in  this  report.  With  the  fre¬ 
quency  response  characteristics  already  determined  by  means  of  the 
digital  computer  procedure  described  in  [2]  and  [3],  all  of  the 
basic  computed  data  and  functions  listed  above  can  be  obtained  in 
a  time  period  of  approximately  30  seconds  on  the  CDC  6600  computer. 
The  computation  time  required  for  determining  the  total  midship 
vertical  bending  moment,  which  is  the  sum  of  the  wave-induced  and 
bow  flare  slam-induced  bending  moments,  is  reduced  significantly 
such  that  the  computation  proceeds  at  a  rate  equivalent  to  170  times 
faster  than  real  time  (based  on  evaluation  corresponding  to  every 
0.05  sec.  of  real  time).  This  is  a  significant  increase  in  computer 
simulation  capability  such  that  rapid  assessment  of  bending  moment 
characteristics,  including  this  type  of  slamming  phenomenon,  can  be 
obtained  simply  and  rapidly  (hence  at  small  cost) .  Thus  a  useful 
todL  for  design  and  analysis  of  ship  structural  response  is  available 
via  computer  simulation  for  this  case,  as  compared  to  model  testing 
and/or  full  scale  tests.  A  flow  chart  illustrating  the  various  pro¬ 
cedures  used  in  computing  time  histories  of  bending  moments  due  to 
wave  action  and  due  to  slamming  is  shown  in  Figure  3. 
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Fig.  3.  Flow  Chart  for  the  Calcu¬ 
lation  of  Slamming  and 
Wave  Bending  Moment  Time 
Histories 
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APPLICATION  TO  BOTTOM  IMPACT  SLAMMING 


The  usual  condition  known  as  slamming  for  ships  in  a  seaway 
is  associated  with  the  emersion  of  the  bow  region  from  the  water 
and  its  subsequent  immersion  upon  re-entering  the  water,  with  an 
associated  impact  force  acting  on  the  ship  during  the  initial  phase 
of  its  re-entry.  Most  studies  of  slamming  phenomenon  have  been 
primarily  concerned  with  the  aspects  of  the  localized  pressure  pulses 
acting  on  portions  of  the  bow  (e.g.  [13]),  which  would  be  responsible 
for  local  damage.  However,  the  concern  in  the  present  investigation 
is  the  resulting  vertical  bending  moment  associated  with  the  struc¬ 
tural  response  of  the  ship,  as  measured  for  reference  purposes  at 
midships.  The  basic  method  of  analysis  will  be  the  modal  technique 
described  previously  in  this  report,  with  the  main  distinction  for 
this  particular  application  being  the  method  of  representation  and 
computation  of  the  exciting  forces  P(x,t)  due  to  ship-wave  interaction. 

The  occurrence  of  slamming  is  associated  with  the  bow  region 
leaving  the  water  and  then  re-entering  at  a  sufficiently  high 
velocity,  which  would  produce  a  high  frequency  local  acceleration 
at  the  bow  as  well  as  local  pressure  fluctuations  and  resulting 
whipping  variations  superimposed  on  the  slowly-varying  direct  wave- 
induced  bending  moment.  The  time  extent  of  the  impact  is  relatively 
short,  and  particular  full  scale  data  such  as  that  in  [14]  for  the 

SS  WOLVERINE  STATE  (a  ship  previously  treated  in  [1],  [2],  and  [7] 

for  wave-induced  bending  moments)  indicates  a  time  duration  extend¬ 
ing  up  to  the  order  of  0.25  sec.  The  basic  mechanism  for  creating 

the  impact  force  in  the  bow  region  is  then  associated  with  the  rapid 

variation  of  effects  that  would  contribute  a  local  vertical  force 
variation  during  re-entry  and  that  may  be  ascribed  to  both  an  added 
buoyancy  as  well  as  an  inertial  force.  The  local  added  buoyancy 
force  can  be  represented  as 


P2(x,t)  =  pgA(zr;x)  (25) 


at  the  local  station  of  interest,  where  A  is  the  cross-sectional 
area  that  varies  in  accordance  with  the  penetration  of  the  section 
into  the  water.  The  inertial  force  is  represented  by 

P1  =  "  5t  (f5wr}  (26) 


where  m  is  the  instantaneous  added  mass  at  the  section  of  interest 
as  the  ship  penetrates  the  water  surface,  wrig  the  relative  velocity 
defined  in  Equation  (19)  ,  and  the  operator  —  is  the  same  as  that 
defined  following  Equation  (18)  .  The  important  point  in  treating 
this  problem  is  to  recognize  that  tae  force  only  occurs  when  the 
ship  section  is  penetrating  the  water  surface  after  re-entering  from 
above,  with  no  force  va]ue  acting  prior  to  the  instant  of  penetra¬ 
tion.  This  reasoning  is  only  applicable  when  considering  the 
determination  of  forces  associated  with  slamming  per  se,  while 
an  assumed  linear  force  variation  is  considered  to  act  on  the 
ship  in  order  to  produce  the  actual  rigid  body  motions,  which  in 
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themselves  are  found  to  be  generally  well  represented  by  linear 
theory  in  spite  of  the  fact  that  a  portion  of  the  ship  may  come 
out  of  the  water. 


Since  the  time  interval  for  generating  the  impact  force 
is  relatively  small  in  a  real  time  sense,  the  time  increment  of 
solution  must  be  reduced  in  order  to  adequately  represent  the 
phenomena  that  take  place.  However,  this  would  tend  to  increase 
computation  time  and  reduce  the  advantage  of  higher  speed  rela¬ 
tive  to  real  time  occurrence,  although  the  possible  extent  of 
reduction  is  not  known  until  much  more  experience  is  obtained 
via  a  large  degree  of  computer  simulation  experimentation.  In 
order  to  carry  out  some  of  the  initial  experiments  on  bottom 
impact  slamming  a  time  increment  (corresponding  to  a  real  time 
solution  interval)  of  0.25  sec.  was  chosen,  and  a  method  establish¬ 
ed  for  predicting  the  possible  occurrence  of  a  slam  by  means  of 
extrapolating  present  position  of  a  ship  section  above  the  instan¬ 
taneous  water  surface  until  a  re-entry  was  expected  (using  the 
value  of  the  relative  vertical  velocity  at  the  present  instant 
of  computed  output) .  The  variation  of  the  local  force  components 
P.  and  P_  defined  by  Equations  (25)  and  (26)  was  determined  during 
tne  time^interval  between  re-entry  and  the  next  time  step  in  the 
computation,  which  would  then  be  a  total  time  interval  of  0.25  sec. 
or  less,  with  a  check  always  being  made  as  to  whether  the  actual 
re-entry  condition  was  achieved  in  order  to  include  the  proper 
force  time  history. 


The  inertial  force  component  was  determined,  in  regard  to 
the  operation  — ,  by  determining  incremental  differences  in  the 
value  of  the  product  mw  and  dividing  by  the  time  interval.  A 
more  detailed  determination  of  the  force  contribution  due  to  the 
term  ^-(inw  )  could  be  obtained  by  means  of  the  following: 
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3  Z 


(27) 


where  an  assumption  of  smaller  time  steps  witKin  the  At  =  0.25  sec. 
interval  is  assumed  in  order  to  allow  for  the  more  accurate  value 
of  r—  that  is  available  from  data  that  has  already  been  tabulated 
(at  i  ft.  intervals) ,  while  assuming  a  relatively  constant  value 
of  w  during  these  intervals  and  also  taking  the  longer  time  step 
(i.e.  0.25  sec.)  in  the  determination  of  wr  .  However,  much  more 
extensive  logic  in  the  computer  program  3t  is  required  in  that  case, 
and  that  approach  was  deferred  to  future  investigations  and  computer 
experiments  of  greater  extent  them  the  present  feasibility  study. 

The  determination  of  time  histories  of  the  local  forces  is 
only  carried  out  for  one  additional  time  step  (0.25  sec.)  beyond 
the  interval  during  which  the  re-entry  of  the  ship  section  occurs. 
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so  that  the  total  force  input  time  history  extent  is  at  most  0 . 5 
sec. ,  which  is  considered  to  be  sufficient  to  represent  the  impact 
time  history  associated  with  this  type  of  slamming  phenomenon. 

While  this  selection  of  time  may  be  somewhat  arbitrary,  and  raises 
the  question  of  possibly  too  large  a  time  extent  of  the  .impact 
force  as  well  as  associated  errors  in  the  resulting  bending  moment, 
a  simple  analysis  illustrates  the  salient  features  of  structural 
response  that  reduces  the  prospects  of  large  errors  due  to  this 
disturbance  time  extent  selection. 

The  modal  analysis  method  represents  the  input  fores  in 
a  weighted  form  given  by  the  function  Q. (t)  defined  by  Equation  (1C) 
with  P(x,t)  *  P.(x,t)  +  P_(x,t).  A  closed  fo.  <r  solution  for  the 
response  q.  (t)  Is  given  by  Equation  (15),  whic’i  can  be  expanded 
into  the  form 
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With  Q. (t)  only  having  a  pulse-like  value  for  a  short  period  of 
time,  Say  t  <  0.5  sec.,  the  integrals  in  Equation  (28)  extend  only 
for  0-*t  ,  and  with  small  system  damping  the  effect  of  the  exponential 
is  negligible.  Thus  the  integral  terms  represent  the  effective 
Fourier  components  of  the  pulse-like  force  Q.  (t)  during  the  time 
extent  t  ,with  the  Fourier  components  being  Chose  associated  with 
the  frequency  w.  (since  X.=w.  for  small  damping  as  3hown  in  Equation 
(16))  and  its  higher  harm&nics  (which  would  be  negligible).  Any 
contributions  from  a  longer  time  extent  for  the  Q  (tj  function  would 
be  expected  to  reflect  the  influence  of  motions  associated  with  the 
lower  frequency  wave-induced  effects  since  the  ship  will  have  pene¬ 
trated  sufficiently  far  into  the  water  that  no  significant  changes 
in  added  mass  and  buoyancy  would  be  occurring  at  a  fast  rate.  Tb-i 
Fourier  components  at  frequency  w.  would  not  be  affected  significantly 
by  such  force  terms  and  hence  no  ippreciable  input  to  the  final 
response  will  occur.  The  response  of  the  system  defined  by 
Equation  (28)  is  then  a  slowly  decaying  transient  at  the  modal 
frequency  w.,  which  is  characteristic  of  the  structural  vibratory 
responses  due  to  slamming,  with  the  response  mainly  in  the  first 
mode,  as  mentioned  previously. 

Computations  of  the  slamming  responses  were  carried  out  for 
the  SS  WOLVERINE  STATE  in  light  load  condition,  as  described  in 
[2] ,  when  proceeding  forward  at  a  speed  of  12  kts.  in  head  seas 
represented  by  a  Sea  State  7  (34  kt.  wind  and  30  ft.  significant 
height) .  The  results  of  motion  computations  to  produce  frequency 
responses  were  obtained  from  [2],  and  this  data  was  used  to  ob¬ 
tain  the  required  kernel  functions  for  use  in  convolution  inte- 


grals  to  obtain  time  histories  of  relative  immersion,  etc.  The 
wave  measurement  reference  point  was  located  35  ft.  ahead  of  the 
ship  FP,  and  the  required  information  to  determine  the  local  forces 
was  obtained  for  the  first  4  stations  of  the  ship  in  the  bow  region 
(stations  1-4),  with  the  ship  assumed  to  be  divided  into  20  stations. 
The  first  mode  frequency  was  selected  as  uk  =  9.42  rad. /sec., 
corresponding  to  the  full  scale  data  reported  in  [14] .  The  mode 
shape  was  taken  to  be  a  simple  parabolic  curve  that  was  similar 
to  that  previously  computed  for  another  Mariner-class  ship,  as 
given  in  [15] .  The  deunping  is  made  up  of  a  structural  damping 
term,  as  well  as  a  term  associated  with  the  spacial  variation  of 
added  mass,  using  the  values  of  the  added  mass  valid  for  high 
frequencies.  This  is  based  upon  the  negligible  influence  of  damp¬ 
ing  due  to  wave  generation  at  the  structural  mode  frequencies,  and 
the  application  of  results  of  ship  motion  strip  theory  (see  [2]  and 
[4]).  Values  of  structural  damping  are  obtained  from  Fig.  2  of 
[4],  so  that  the  damping  parameter  C.  in  Equation  (6)  is  defined 
by  1 
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where  A',,  is  the  local  sectional  vertical  added  mass  (determined 
in  this  case  for  the  high  frequency  limit) .  The  value  of  total 
damning  for  the  present  case  (first  model  response)  was  found  to 
be  _1  =  0.13,  with  the  structural  damping  portion  corresponding  to 

o.osi. 

An  analysis  of  the  wave-induced  midship  bending  moment  for 
the  WOLVERINE  STATE  in  the  frequency  domain  (by  using  the  methods 
and  results  of  [2])  yields  the  response  amplitude  operator  (ampli¬ 
tude  of  bending  moment  per  unit  wave  amplitude)  given  in  Figure  4 
for  this  case.  The  power  spectrum  of  this  bending  moment  resulting 
from  assuming  a  wave  spectrum  corresponding  to  Sea  State  7  is 
shown  in  Figure  5  (with  total  spectral  area  =  o2,  where  o  is 
the  rms  value) ,  with  a  resultant  rms  bending  moment  value  of 
3.55  x  104  ft. -tons.  The  amplitude  response  in  Figure  4,  together 
with  phase  information,  was  used  to  obtain  the  kernel  function  for 
the  wave-induced  bending  moment.  This  allows  generation  of  time 
domain  records  of  this  bending  moment  component  for  a  given  wave 
record  input  that  could  be  compared  with  the  slam-induced  bending 
moment  time  histories  via  the  computer  simulation  technique  de¬ 
scribed  above. 

The  results  of  the  computations  of  slamming  responses  were 
obtained  for  an  extensive  run  in  time,  where  the  occurrence  of 
slamming  was  indicated  directly  in  the  computer  output.  As  men¬ 
tioned  previously  all  computations  were  carried  out  during  an 
equivalent  time  interval,  corresponding  to  real  time  of  0.25  sec., 
while  the  solution  of  the  modal  responses  from  Equation  (6)  was 
carried  out  by  means  of  the  digital  technique  described  by  Equa¬ 
tion  (17)  with  a  sampling  time  of  every  0.05  sec.  This  was  done 


m%t  r«4./t#c. 

Fig.  4.  Midship  Wave-Induced  Bending 
Moment  Response  Amplitude 
Operator,  SS  WOLVERINE  STATE, 
V  =  12  kts..  Head  Seas 
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Fig.  5.  Midship  Wave-Induced  Bending  Mo¬ 
ment  Power  Spectrum,  SS  WOLVERINE 
STATE,  V  =  12  kts..  Sea  State  7 
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Fig.  6.  Time  Histories  of 
Bending  Moments  and 
Associated  Forces, 
WOLVERINE  STATE,  V  = 
12.0  kts..  Sea  State  7 
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in  order  to  properly  obtain  the  high  frequency  responses  of  the 
system,  and  it  was  also  the  basis  for  determining  the  responses  due 
to  bow-flare  slamming  treated  in  the  previous  section. 

A  typical  output  is  shown  in  Figure  6,  which  provides  a 
time  history  of  the  wave  (as  measured  35  ft.  ahead  of  the  ship)  , 
the  wave  induced  bending  moment,  and  the  bending  moment  due  to 
slamming.  Since  the  system  has  relatively  small  damping  the  re¬ 
sults  of  earlier  slams  are  not  sufficiently  decayed,  and  there  is 
often  a  mixture  of  contributions  to  the  total  slam  output.  The 
data  in  Figure  6  also  shows  the  separate  contributions  of  the  two 
force  terms  at  station  2  which  contribute  toward  the  slam  response 
in  this  case.  The  results  obtained  for  the  bending  moments  associated 
with  the  bottom  impact  slam  are  much  too  large,  since  they  are  the 
same  order  of  magnitude  as  the  wave-induced  bending  moments  and  that 
is  not  the  experience  generally  indicated  in  the  results  of  [14] . 
However  the  relative  vertical  velocities  in  the  bow  region  associated 
with  the  occurrence  of  slamming  were  found  to  be  in  the  range  of 
12-20  ft. /sec.  ,  which  is  in  general  agreement  with  such  results  in 
[14]. 


The  large  responses  can  be  easily  ascribed  to  the  effect  of 
too  large  an  input  force,  mainly  due  to  the  inertial  term  P,  given 
in  Equation  (26) .  This  quantity  is  affected  by  the  value  of  At 
chosen  for  the  determination  of  the  time  derivative,  and  in  addition 
the  force  magnitude  is  based  upon  the  very  simple  model  indicated 
by  Equations  (25)  and  (26) .  It  is  known  that  other  effects  are 
present  which  would  account  for  the  influence  of  the  actual  deadrise 
form  of  the  section,  which  would  also  account  for  local  water  rise 
due  to  the  impact,  and  might  possibly  require  consideration  of  the 
effect  of  "cushioning"  due  to  entrapped  air,  etc.  All  of  these 
effects  require  a  more  sophisticated  model  for  representing  the 
forces,  which  goes  beyond  even  the  more  accurate  representation  for 
the  simplified  inertial  reaction  shown  in  Equation  (27) .  It 
does  not  appear  to  be  difficult  to  carry  out  the  determination  of 
the  bending  moment  response  via  the  present  methods  as  long  as 
an  accurate  force  input  is  included  within  the  overall  representa¬ 
tion.  The  time  scale  for  the  present  simulation  was  about  80  times 
faster  than  real  time;  which  reflects  the  increased  computation 
time  required  by  reducing  the  time  increment  from  0.50  sec.  to  0.25 
sec.  in  this  phase  as  compared  to  the  case  of  bow-flare  slamming 
treated  earlier.  Since  the  present  program  was  a  feasibility  study 
for  computer  simulation,  and  complicated  force  mechanisms  were  not 
to  be  covered  in  detail,  the  basic  capability  of  this  simulation 
is  indicated  by  the  present  results.  However,  more  detailed  analysis 
and  representation  of  the  local  force  variation,  which  can  be  ob¬ 
tained  from  many  sources  and  which  must  then  be  converted  into  a 
useful  computational  form  for  evaluation  of  the  forces  within  a 
small  time  increment,  is  the  main  task  required  for  more  precise 
slam-induced  bending  moment  evaluations.  This  lies  beyond  the 
present  investigation  requirements  and  remains  as  a  future  task 
in  continuation  of  ship  structural  response  evaluation  via  computer 
simulation.  However  the  basic  computational  procedures  for  evaluat¬ 
ing  the  slam-induced  bending  moments  are  essentially  the  same  as 
for  the  bow  flare  case,  and  are  thus  outlined  in  a  flow  chart  form 
very  similar  to  that  in  Figure  3. 
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APPLICATION  TO  SPRINGING 


When  considering  the  case  of  springing,  the  assumption  is 
made  that  there  is  negligible  influence  of  the  ship  rigid  body 
motions  on  this  particular  phenomenon  and  that  the  external  force 
acting  on  the  vibratory  system  is  only  that  due  to  the  waves  in 
their  interaction  with  the  ship.  The  local  wave  force  acting  on 
a  section  of  the  ship,  which  represents  the  force  P(x,t)  in  the 
model  equation  model,  is  given  by 
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is  the  surface  wave  elevation  and 
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The  quantity  N'  is  the  local  sectional  vertical  damping  force 
coefficient,  alfid  E  is  the  mean  section  draft  that  is  approximated 
by 
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where  H  is  the  local  section  draft  and  C  is  the  local  section  area 
coefficient.  This  expression  is  obtained  from  [2]  and  represents 
the  result  of  the  usual  technique  of  strip  theory  to  express  the 
vertical  wave  force  acting  on  a  ship. 

In  determining  the  quantity  Q. (t)  defined  in  Equation  (10) 
it  is  necessary  to  weight  the  wave  forces  in  terms  of  the  mode 
shape  and  then  integrate  the  result  over  the  ship  length.  The 
computations  must  be  carried  out  over  a  range  of  wave  lengths  that 
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would  extend  down  to  small  values  relative  to  the  ship  length  in 
order  to  excite  the  higher  frequencies  that  correspond  to  the  ship 
structural  modes,  with  the  main  influence  being  that  of  the  first 
mode  of  vibration,  as  discussed  previously.  The  particular  require¬ 
ments  for  carrying  out  such  an  integration  in  an  accurate  fashion 
by  numerical  means,  as  well  as  the  results  of  application  to  a 
specific  ship  case,  will  be  discussed  in  a  later  portion  of  this 
section  of  the  report. 

In  order  to  obtain  time  histories  for  the  case  of  springing, 
it  would  be  necessary  to  obtain  values  of  the  wave  time  history 
and  its  next  two  time  derivatives  at  various  stations  along  the 
hull,  with  proper  representation  of  the  phase  differences  due  to 
location.  All  of  these  terms  must  be  related  to  the  initial  single 
wave  motion  time  history  as  measured  at  a  reference  point  ahead  of 
the  ship,  and  this  would  require  at  least  one  kernel  function  for 
eac'i  station  (assuming  20  stations  for  the  ship)  as  well  as  the 
requirements  of  a  highly  oscillatory  kernel  function  for  represent¬ 
ing  only  a  phase  shift  frequency  response (see  [16]).  When  coupled 
with  the  different  spacial  variations  due  to  differing  geometry 
at  each  of  the  stations  along  the  ship,  the  computational  complexity 
associated  with  this  procedure  is  very  evident.  In  order  to  pro¬ 
ceed  with  a  time  domain  representation  of  the  springing  response, 
it  is  then  necessary  to  directly  obtain  frequency  responses  for 
the  bending  moment  associated  with  the  springing,  which  is  obtained 
from  the  frequency  response  representation  of  the  basic  modal 
response,  i.e.  Equation  (6).  This  will  produce  a  resulting  kernel 
function  for  the  bending  moment  due  to  springing,  which  is  expected 
to  reflect  the  sharply  tuned  vibratory  response  associated  with  a 
lightly  damped  second  order  differential  equation. 

The  solution  for  the  vibratory  response  in  the  frequency 
domain  is  obtained  by  representing  the  forcing  function  Q. (t)  in 
the  form  1 


Q^t)  -  Qoej“et 


(35) 


so  that  the  steady  state  solution  of  Equation  (6)  is  given  by 

0oe)“et 


Ki”“e2wi+^weCi 


(36) 


in  complex  form,  from  which  the  amplitude  and  phase  can  be  obtained 
as  a  function  of  the  frequency  u>  .  It  can  be  easily  shown  that 
this  solution  leads  to  the  same  analytical  result  for  the  vertical 
elastic  deflection  as  given  in  [4] .  The  midship  bending  moment  is 
given  by  the  present  analysis  as 
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where  the  significant  response  is  given  only  for  i«l,  correspond¬ 
ing  to  the  first  mode  of  vibration  of  the  ship  structure,  and  this 
form  can  be  expressed  in  terms  of  an  amplitude  and  phase  relative 
to  a  wave  time  history  reference  at  a  particular  location  with 
respect  to  the  ship. 

An  application  of  this  approach  was  made  for  a  ship  selected 
to  represent  a  200,000  dwt.  tanker  which  was  based  on  the  Series  60, 
block  coefficient  0.80,  form  that  was  previously  analyzed  for  wave- 
induced  bending  moments  in  [2] .  The  ship  had  a  displacement  of 
250,000  tons,  with  a  length  of  1,100  ft.  and  other  geometric  para¬ 
meters  corresponding  to  that  particular  Series  60  ship  form.  The 
first  mode  natural  frequency  was  selected  as  w.  =  3.0  rad. /sec. , 
and  the  mode  shape  was  taken  to  correspond  to  the  same  basic  para¬ 
bolic  form  as  used  for  the  WOLVERINE  STATE  (obtained  from  [15])  but 
properly  scaled  to  correspond  to  the  dimensions  of  the  ship  selected. 
The  forward  speed  was  assumed  to  correspond  to  Froude  No.  =  0.15, 
which  was  a  speed  of  16.7  knots. 

The  wave  force  defined  by  Equation  (30)  was  determined  for 
a  series  of  waves  corresponding  to  the  range  0 .  06<id^2. 265 ,  thereby 
covering  the  normal  rigid  body  ship  response  region  as  well  as  the 
range  of  short  waves  that  would  have  encounter  frequencies  close 
to  first  vibration  mode  frequency.  The  computations  of  P(x,t)  were 
carried  out  by  the  use  of  the  digital  computer  program  in  [3]  ,  using 
the  complete  frequency-dependent  coefficient  values  determined  from 
that  program  for  each  frequency  (at  higher  frequencies  the  expected 
limits  of  the  various  coefficient  terms  such  as  N'  and  A'  were 
properly  obtained)  .  The  integration  over  the  shi£  length,  with  the 
mode  shape  weighting  X.(x),  to  produce  the  function  Q. (t)  defined 
in  Equation  (10)  was  carried  out  using  Filon's  methodiof  integration 
[17]  in  order  to  properly  account  for  the  influence  of  the  short 
waves  relative  to  the  station  spacing  corresponding  to  a  20  station 
representation . 

The  analysis  described  above  for  the  case  of  springing  is 
based  on  the  modal  response  model,  with  the  resulting  bending  moment 
due  to  springing  represented  by  Equation  (37).  This  result  is  ob¬ 
tained  for  the  general  case  where  it  is  assumed  that  responses  are 
characterized  directly  in  terms  of  a  sum  of  separate  results  at  each 
of  the  modal  frequencies,  as  a  result  of  transient  impulsive  inputs, 
and  that  form  is  also  inherent  in  the  results  of  the  analysis  in  [5] . 

A  simple  examination  of  the  consequence  of  the  present  analysis, 
as  well  as  those  exhibited  in  [5]  ,  shows  a  significant  peak  in  the 
springing  bending  moment  frequency  response  in  the  region  where  verti¬ 
cal  ship  motions  are  usually  present.  Since  such  motions  have  been 
neglected  in  the  present  case,  there  is  a  question  as  to  the 
significance  and  validity  of  any  large  response  indicated  at 
those  frequencies  due  to  the  springing  phenomenon,  which  would 
not  be  expected  to  excite  structural  modes  at  the  lower  frequencies 
where  wave-induced  bending  moments  and  ship  motions  are  predominant. 

An  examination  of  the  results  for  the  problem  of  springing 
given  in  [4]  shows  a  somewhat  different  behavior  in  the  low  frequency 
range,  since  the  inertial  reactions  associated  with  the  ship  struc¬ 
tural  deflection  accelerations  introduce  terms  proportional  to  u  2 
which  significantly  reduce  any  possible  contribution  at  the  low 
frequencies  (while  still  neglecting  the  ship  rigid  body  motions) . 
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While  the  influence  of  such  an  inertial  reaction  at  frequencies 
higher  than  the  first  mode  resonance  would  tend  to  amplify  the 
bending  moment  response  for  this  mode  somewhat,  the  oscillatory 
decay  at  higher  frequencies  of  the  wave  excitation  forces  Q,(t)  in 
the  deflection  representation  given  by  Equation  (36)  will  reduce 
the  springing  bending  moment  response  significantly  in  this  higher 
frequency  range.  Thus  the  form  of  the  frequency  response  for  the 
midship  bending  moment  due  to  springing  should  be  more  properly 
represented  by 
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with  a  rapid  asymptotic  decay  in  frequency  beyond  the  first  mode 
value.  The  precise  form  of  the  wave  excitation  forces  represented 
by  P(x,t)  and  the  resulting  Q.  (t)  is  not  known,  although  the  value 
given  by  Equation  (30)  is  usea  here  and  has  also  been  applied  in  [4] 
and  [5]  .  It  is  known  that  a  rapid  decay  of  these  quantities  with 
frequency  will  occur,  so  that  the  frequency  response  form  beyond 
the  first  mode  resonance  will  not  be  significantly  altered.  A  more 
complete  discussion  and  analysis  of  wave  forces  for  short  wavelengths 
relative  to  the  ship  length  will  be  given  in  the  next  section  of 
this  report. 

The  frequency  response  functions  were  determined  from  Equations 
(36)cand  (38),  where  the  damping  value  for  this  case  corresponded 
to  _1  =  0.044,  with  the  amount  0.008  due  to  structural  damping. 

Figures  7  and  8  represent  the  midship  bending  moment  frequency 
responses,  from  which  a  time  domain  kernel  function  corresponding 
to  the  wave  reference  point  located  30  ft.  ahead  of  the  FP  was  cal¬ 
culated.  This  springing  bending  moment  kernel  function  is  shown  in 
Figure  9,  and  the  kernel  function  for  the  wave-induced  midship  bend¬ 
ing  moment  is  shown  in  Figure  10. 

Computations  of  time  history  responses  of  these  bending  moments 
were  carried  out  for  this  ship  moving  at  16.7  knots  in  a  Sea  State  5 
which  corresponds  to  a  22  kt.  wind  speed  and  a  significant  height 
of  10  ft.  The  filter  representing  the  encountered  wave  motion 
spectrum  was  established  by  the  approximation  technique  described 
in  [12] ,  which  characterizes  the  major  properties  and  form  of  the 
spectrum,  although  there  may  not  be  exact  correspondence  in  the 
spectral  ordinates  at  all  frequencies,  especially  the  higher 
values  that  would  be  of  interest  in  the  present  problem.  The 
wave  spectral  ordinate  at  the  frequency  u>  =  3.0  is  0.8  ft.2- 
sec.  from  the  filter,  while  the  theoreticll  value  for  an  idealized 
wave  spectrum  is  0.5  ft.2-sec.,  so  that  a  large  springing  response 
will  be  obtained  here  in  the  present  simulation  as  compared  to  a 
value  predicted  for  an  idealized  wave  spectrum. 

The  computed  time  history  responses  of  the  wave  induced 
bending  moment  and  the  springing  bending  moment  are  shown  in 
Figure  11  together  with  the  wave  record  corresponding  to  this 
Sea  State  5.  A  time  increment  of  .25  sec.  (real-time)  was  chosen 
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Fig.  8.  Midship  Wave-Induced  Bending 
Moment  Response  Amplitude 
Operator,  200,000  dwt.  Tanker, 
V  =  16.7  kts. 


Fig.  7.  Midship  Springing  Bending  Moment 
Response  Amplitude  Operator, 
200,000  dwt.  Tanker,  V  ■  16.7  kts. 
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Fig.  9.  Kernel  Function  for  Midship 
Springing  Bending  Moment, 

200,000  dwt.  Tanker,  V  ■  16.7  kts. 


m 


m 


Fig.  JO.  Kernel  Function  for  Wave-Induced  Midship 
Bending  Moments,  200,000  dwt.  Tanker, 

V  =  16.7  kts 
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Fig.  11.  Midship  Bending  Moment  Time  Histories, 

200,000  dwt.  Tanker,  V  =  16.7  kts. 

for  this  phase  of  the  study  to  properly  handle  the  higher  fre¬ 
quencies  present  in  the  springing  phenomena.  The  time  histories 
shown  in  Figure  11  clearly  demonstrate  the  expected  sharply  tuned 
vibratory  response  associated  with  springing.  The  correctness  of 
the  magnitude  of  these  moments  can  be  ascertained  by  a  comparison 
of  the  rms  values  obtained  by  an  analysis  of  the  time  histories 
over  a  sufficiently  long  period  of  time  with  the  rms  values  ob¬ 
tained  from  the  areas  under  their  respective  power  spectra.  The 
power  spectra  used  in  this  comparison  is  based  on  the  wave  generat¬ 
ing  digital  filter  rather  than  on  an  idealized  wave  spectra,  since 
the  purpose  of  the  comparison  is  to  validate  the  feasibility  of 
using  this  technique  of  simulating  the  time  history  wave  induced 
bending  moment  and  springing  bending  moment  experienced  by  ship 
in  a  random  sea.  The  rms  values  of  the  wave-induced  bending  moment 
are  1.3  x  105  ft. -tons  and  1.5  x  105ft.-tons  from  the  time  history 
analysis  and  frequency  domain  analysis,  respectively.  The  corres¬ 
ponding  rms  values  for  the  springing  bending  moment  are  4.4  x  10 5 
ft. -tons  and  4.1  x  10s  ft. -tons.  This  close  agreement  demonstrates 
the  feasibility  of  using  the  time  domain  representation  of  the 
springing  response.  The  fact  that  the  springing  response  has  the 


larger  nns  value  is  solely  the  result  of  the  fall*-off  characteristics 
of  the  particular  digital  filter  chosen  for  the  study.  A  more 
exact  digital  filter  representing  the  true  wave  spectra  is  easily 
obtainable  if  extreme  resolution  in  duplicating  frequency  character¬ 
istics  is  required  in  this  range. 

Another  feature  of  this  tiirs  domain  analysis  is  the  fact 
that  only  two  convolution  integrals  are  required  for  the  sirulation. 
As  a  result,  the  computation  of  response  time  histories  can  be 
run  150  times  faster  than  real  time.  The  flow  chart  representing 
the  computational  procedures  used  in  this  time  history  simulation 
for  bending  moments  due  to  springing  and  wave-induced  (low  fre¬ 
quency)  effects  is  given  in  Figure  12. 


Fig.  12.  Flow  Chart  for  the  Calculation 
of  Springing  and  Wave  Bending 
Moment  Time  Histories 

WAVE  FORCE  AT  HIGH  FREQUENCY 

As  mentioned  previously,  as  well  as  in  the  work  of  [4] 
and  (5],  the  dependence  of  springing  response  on  the  wave  force 
distribution  along  the  ship  hull  in  the  high  frequency  range  where 
this  phenomenon  occurs  requires  an  examination  of  these  forces 
in  that  range.  The  force  expression  used  is  valid  for  wavelengths 
that  are  of  the  order  of  the  ship  length  (kL=0(l)),  where 

k=  -  —  )  and  hence  large  compared  to  the  cross  sectional 
dimensionS  (see  (18])  and  [19]),  while  th^  important  wavelengths 
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in  the  case  of  springing  are  of  the  order  of  the  ship  beam  or 
draft. 

In  order  to  obtain  further  insight  into  the  possible  varia¬ 
tion  of  the  force  at  high  frequencies,  a  limited  analysis  was 
carried  out  to  derive  the  wave  force  under  the  initial  assumption 
of  high  frequency.  The  analysis  is  restricted  to  the  case  of  zero 
forward  speed,  for  simplicity,  and  the  assumption  of  treating  a 
two-dimensional  problem  in  the  cross-flow  plane  for  head  seas  is 
made.  This  latter  assumption  is  based  upon  the  form  of  the  wave 
potejjij^al  and  associated  interaction  potential  being  proportional 
to  ex  along  the  hull,  which  has  already  been  shown  to  be  appli¬ 
cable  for  submerged  bodies  in  short  waves  (kR=0(l),  with  R  the 
radius).  The  problem  is  established  in  the  two-dimensional  cross- 
flow  plane  for  a  simple  body,  viz.  a  circle,  as  shown  in  the 
sketch  below 


z 'plane 


The  incident  wave  potential  is  represented  by 

♦„  -  «celt*,ei<,tx-“t,  (39) 

where  the  x-coordinate  is  along  the  body  length  (to  be  distinguished 
from  the  x' -coordinate  in  the  body  sketch  above) .  The  two-dimen¬ 
sional  Laplace  equation  is  assumed  to  apply  for  the  cross-flow 
potential  and  this  is  solved  by  initially  mapping  the  flow  region 
into  another  plane,  i.e.  by  use  of  the  mapping  transformation 


0;  y'  ■  0,  outside  the  body 


for  the  high  frequency  case.  The  vertical  velocity  on  the  cir¬ 
cular  boundary  is  known  from  the  wave  potential, 

v  -  -  akc  ei(kx-u,t)ekRsin  0  I 


using  the  polar  coordinate  form  in  the  z' -plane.  Thus  a  mixed 
boundary  value  problem  is  established  for  the  disturbance  potential 
and  this  is  solved  by  initially  mapping  the  flow  region  into 
another  plane,  i.e.  by  use  of  the  mapping  transformation 


(42) 
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The  velocities  in  the  two  planes  are  related  by 

-  u'+iv'  ■  (-u+iv)  |"lt  - - -  I 

*-  /  c^-R^  ■* 


(43) 


where  the  primed  (')  quantities  are  those  in  the  c -plane,  and  the 
resulting  boundary  value  problem  for  the  Laplace  equation  in  this 
plane  is  shown  in  the  sketch  below 


C-plane 


$■0 
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On  the  plate  5,  where  |(|<R,  so  that  v'=v,  which  is  expressed  by 


V  =  akcei(kx‘u,t)e"kR/l‘U/R)T 


(44) 


The  mixed  boundary  value  problem  can  be  converted  into  a 
singular  integral  equation  by  applying  Green's  theorem,  leading  to 
the  Poisson  formula  for  the  half-plane,  and  that  in  turn  by  dif¬ 
ferentiation  and  integration  operations  with  respect  to  the  vari¬ 
ables  n '  and  £  can  be  shown  to  lead  to  the  equation 
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where  d  is  a  constant.  The  solution  of  this  equation, 
it  is  bounded  at  both  ends  (-R  and  R)  ,  is  given  by 
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which  is  expressed  by 
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(see  [20]  for  discussion  of  solutions  of  singular  integral  equa¬ 
tions  ) . 

The  vertical  force  acting  on  the  section  is  given  by 


-pR  8^n  0  dO  +  p  ft^ 


where  the  part  due  to  $  represents  the  direct  incident  wave  force, 
which  is  primarily  hydrostatic,  while  the  part  due  to  $  is  the  effect 
of  wave-body  interaction.  Using  the  previous  expressions  and  carry¬ 
ing  out  the  indicated  integration  of  the  pressure  represented  by 
the  operations  in  Equation  (49) ,  leads  to 


F1  -  -  pR 


sin  e  de 


-  ipgaRei(kx-“t) 
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The  integral  can  be  expressed  as 
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which  can  be  evaluated  as 
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given  by**  0thei  f°rCe  component'  due  to  wave-body  interaction,  is 


p2  =  P  d5 

■  -H 

which  can  be  reduced  to  a  single  integral  expression  given  by 
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The  total  wave  force  equal  to  the  sum  of  and  F2  is  then  given  by 


F/pgaR=  2“"  (kR)  Po(kR)"Lo(kR)l 


which  will  be  valid  for  high  frequencies  according  to  the  previous 
derivation.  Another  approach  was  also  used  to  treat  this  particular 
problem  based  on  establishing  an  analysis  that  treats  the  high 
frequency  problem  by  the  method  of  matched  asymptotic  expansions 
[22] ,  and  the  same  exact  result  was  obtained  [23] .  However  the  de¬ 
tailed  consideration  of  the  effects  of  high  frequency  led  to  an 
examination  of  the  phenomena  associated  with  the  flow  in  the  region 
of  the  "corners"  representing  the  intersection  at  the  ends  of  the 
section  with  the  free  surface.  The  corner  problem  must  be  solved, 
subject  to  a  complicated  differential  equation  that  represents  the 
flow  conditions  there,  and  that  difficult  problem  must  be  treated 
in  order  to  obtain  a  complete  solution  for  the  present  case  of 
high  frequency  waves.  Since  that  is  beyond  the  scope  of  the 
present  investigation,  no  further  consideration  was  given  to  that 
problem. 

In  order  to  compare  the  results  of  the  present  analysis, 
Equation  (54) ,  with  the  conventional  value  obtained  from  applying 
strip  theory  to  this  semicircle  case,  the  results  of  Equation  (30) 
were  found  for  this  configuration.  The  hydrostatic  and  added  mass 


Fig.  13.  Comparison  of  Theoretical  Ex¬ 
pressions  for  Sectional  Wave 
Force  on  a  Semicir  le 

section  form  and  the  associated  effect  of  the  sinusoidal  functions 
of  kx  that  are  required,  will  not  result  in  any  significant  force 
differences.  The  variation  in  force  amplitude  for  different,  wave¬ 
lengths  will  be  mainly  manifested  by  the  influence  of  the  e1  var¬ 
iation  within  the  integrations,  with  only  a  small  effect  of  the 
other  weighting  functions.  This  may  possibly  be  altered  as  a  result 
of  including  forward  speed  effects,  which  would  involve  more  analytic¬ 
al  complexity,  or  by  the  influence  of  the  corner  solution  discussed 
previously.  However  no  further  consideration  can  be  given  to  this 
problem  at  this  time  since  it  is  beyond  the  scope  of  the  present 
investigation,  and  it  remains  as  a  task  for  future  detailed  studies. 

At  the  present  time  though  it  appears  that  from  the  known  available 
results,  as  shown  above,  that  the  strip  theory  values  for  wave 
force  can  still  be  used  in  the  analysis  of  springing  until  such  time 
as  a  more  comprehensive  theory  (or  results  of  experiments  on  high 
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frequency  wave  force  measurements)  shows  significant  differences 
from  those  values. 


DISCUSSION  AND  CONCLUDING  REMARKS 

The  results  obtained  in  the  preceding  sections  show  that 
computer  simulation  techniques  can  be  applied  to  determine  the 
vibratory  structural  responses  of  ships  in  waves  at  a  computational 
rate  significantly  faster  than  that  required  in  the  equivalent 
real  time  of  the  actual  ship  mocion.  The  resulting  time  domain 
outputs  can  be  displayed  for  detailed  examination  together  with 
information  on  local  force  values,  the  time  history  of  the  refer¬ 
ence  waves  that  are  basically  responsible  for  bringing  about  the 
ship  response,  phase  relations,  etc.  In  addition  the  time  domain 
outputs  can  be  continuously  processed  on-line,  i.e.  simultaneously 
while  they  are  being  generated,  in  order  to  provide  statistical 
measures  such  as  rms  values,  peak  values,  level  exceedance  probabil¬ 
ity,  and  other  desired  properties  of  a  time  history  record  output 
just  as  is  presently  being  done  with  full  scale  records  from  ships 
at  sea  (e.g.  [14]),  but  with  the  added  advantage  of  obtaining  the 
resists  at  a  much  faster  rate  than  real-time.  All  of  these  results 
are  obtained  using  a  large  high  speed  digital  computer  with  programs 
written  in  the  FORTRAN  language.  Thus  a  general  utility  can  be 
achieved  in  determining  this  type  of  structural  response  by  applica¬ 
tion  of  commercially  available  general  purpose  computers  rather 
than  being  restricted  to  a  particular  special  purpose  computer  or 
class  of  computers  (such  as  hybrid  systems,  etc.  that  were  initially 
applied  in  [1]). 

The  capability  of  achieving  these  results,  with  values  that 
appear  to  be  close  to  those  obtained  in  experimental  studies,  is 
well  demonstrated  in  the  case  of  bow-flare  slamming,  where  the 
physical  model  of  the  force  mechanism  is  fairly  well  established 
and  it  is  also  a  relatively  continuous  type  of  force  variation.  In 
the  case  of  bottom  impact  slamming  however,  the  force  mechanism  used 
herein  was  a  simple  model  that  also  was  limited  in  its  simulation 
due  to  the  time  step  restrictions  imposed  in  the  computations.  A 
more  realistic  representation  of  the  computations,  such  as  that 
given  in  Equation  (27) ,  would  appear  to  provide  a  better  force 
expression  that  could  reduce  the  influence  of  the  increment  of  At 
used  in  the  determination  of  time  rate  change.  However  even  the 
basic  force  mechanism  used  here  would  not  be  expected  to  be  fully 
appropriate,  since  a  relatively  flat  bottom  would  produce  very 
large  changes  in  the  time  derivative  that  would  not  be  realistic. 
Thus  the  task  required  In  the  case  of '’simulating  bottom  impact 
slamming  responses  is  a  better  representation  of  the  local  force 
associated  with  the  entry  of  a  ship  section  into  the  water.  There 
are  many  possible  models  to  use  for  this  phenomenon,  and  they  must 
be  evaluated  for  their  applicability  since  accurate  representations 
of  certain  models  of  the  flow  field,  such  as  the  inclusion  of 
entrapped  air  and  its  effect  on  the  local  water  surface  deformation, 
could  not  be  reasonably  expected  to  be  proper  (or  even  capable  of 
adequate  mathematical  representation)  when  considering  the  real 
case  of  relative  motion  involving  irregular  wave  forms,  possible 
spray  and  local  air  content  in  the  sea  water,  etc.  While  previous 
studies  of  impact  on  the  water  surface,  when  related  to  slamming, 


have  been  mainly  concerned  with  determining  local  pressure  time 
histories  at  particular  points  (such  as  along  the  centerline,  etc.) 
or  some  representation  of  spacial  pressure  distributions  (with 
knowledge  of  their  time  history  variations),  the  present  problem 
requires  determining  the  vertical  force  time  history.  The  pre¬ 
sence  of  the  "smoothing"  effect  of  spacial  integration,  which  might 
be  accomplished  by  analytic  means  initially,  will  tend  to  reduce 
the  complexity  of  the  computations  as  well  as  provide  a  more 
accurate  total  force  representation. 

The  present  computer  experiments  have  demonstrated  the 
feasibility  of  obtaining  structural  responses  by  the  methods 
developed  herein,  and  it  then  remains  to  apply  more  precise  local 
force  terms  as  an  input.  The  most  promising  physical  and  mathem¬ 
atical  models  for  such  forces  must  be  analyzed  to  determine  their 
applicability  to  this  present  problem,  which  must  include  evaluating 
the  computational  efforts  necessary  for  the  force  representation 
vis-a-vis  the  extent  of  final  structural  response  output  relative 
to  experimental  results,  in  order  not  to  compromise  any  of  the 
benefits  possible  in  producing  outputs  at  a  rate  significantly  faster 
than  real-time.  This  task  is  recommended  as  an  important  element 
to  be  accomplished  in  any  further  investigations  of  ship  computer 
response  studies. 

For  the  case  of  springing,  the  only  questionable  aspect 
appears  to  be  the  accurate  representation  of  the  wave  force  along 
the  hull  for  short  wavelengths.  The  analytical  studies  described 
in  the  preceding  section  of  this  report  did  not  seem  to  indicate 
any  significant  alteration  in  the  force  value  relative  to  the  value 
obtained  by  usual  strip  theory  methods.  However,  the  effect  of 
forward  speed  should  be  examined,  and  similarly  any  possible  in¬ 
fluence  of  the  "corner  problem"  described  in  the  last  section 
should  also  be  studied,  in  order  to  cover  the  possible  variations 
in  the  hydrodynamic  wave  force  that  would  influence  the  magnitude 
of  springing  response.  Incorporation  of  such  results  into  the 
available  mathematical  and  computational  model  can  be  easily 
accomplished,  resulting  in  a  useful  tool  for  predicting  the  mag¬ 
nitude  and  properties  of  that  type  of  structural  response. 

The  form  of  representation  of  springing  response  can  be 
given  in  terms  of  spectral  output,  in  the  frequency  domain,  or 
alternatively  time  histories  can  be  provided.  Since  the  fundamental 
phenomena  are  linear,  the  output  in  the  form  of  spectra  and  asso¬ 
ciated  statistical  properties  such  as  rms ,  significant  amplitude, 
etc.  are  more  easily  determined  in  that  manner.  Time  history  out¬ 
puts  are  also  useful  for  correlation  purposes,  examining  the  decay 
properties  of  the  vibrations,  and  other  detailed  studies  of  a 
representative  record  of  response.  Another  useful  application  of 
time  history  outputs  is  for  comparison  with  model  experimental 
data,  since  magnitude  and  phase  for  regular  waves  (and  precise 
time  history  duplication)  can  be  easily  compared  to  determine  the 
validity  of  the  theory  in  detail. 

Some  consideration  should  also  be  given  to  conducting  verifi¬ 
cation  model  experiments  for  the  vibratory  structural  responses, 
using  a  segmented  model  that  contains  sufficient  flexibility  to 
provide  vibratory  structural  response  data.  The  model  should  have 
all  of  its  structural  and  inertial  properties  fully  defined  for 


purposes  of  direct  comparison  to  theoretical  computed  results.  An 
interchangeable  bow  region  can  be  built  so  that  it  can  reflect  the 
form,  flare,  bottom  contour,  etc.  that  would  be  appropriate  to  such 
ships  as  container  vessels,  normal  cargo  ships,  tankers,  etc.  Thus 
different  tests  can  be  carried  out  in  waves  to  represent  the  effect 
of  slamming  (bow  flare  and  bottom  impact)  when  the  ship  has  large 
motions  in  the  waves,  as  well  as  for  springing  by  testing  in  short 
waves.  Careful  design  of  the  model  and  test  programs  will  allow 
covering  these  different  and  important  structural  response  phenomena, 
and  wxll  provide  a  direct  and  useful  validation  of  results  obtained 
from  computer  simulation  studies. 

The  various  analytical  and  computational  procedures  de¬ 
scribed  above  can  be  readily  accomplished,  and  represent  a  logical 
continuation  of  the  efforts  originally  established  for  the  ship 
computer  response  project.  The  final  product  will  be  a  digital  com¬ 
puter  program,  similar  to  that  described  in  [2]  and  [3],  that  will 
enable  users  to  obtain  outputs  for  determining  structural  responses 
of  vibratory  nature  such  as  those  resulting  from  slamming  and  spring¬ 
ing.  Thus  a  complete  evaluation  of  the  major  structural  loads  and 
responses  of  a  ship  in  a  seaway  can  then  be  predicted  in  a  relatively 
short  time  by  use  of  that  type  of  program,  in  conjunction  with  that 
in  [3]  ,  as  an  important  tool  in  design  and  analysis  studies.  De¬ 
pending  upon  the  extent  of  validation  desired,  the  experimental  pro¬ 
gram  described  above  can  also  be  carried  out  as  an  adjunct  to  the 
computer  program  development  in  order  to  assess  the  reliability  of 
the  analytical/computational  procedures  in  providing  realistic  values 
of  ship  structural  loads. 
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